The master equation for the reduced open-system dynamics, 
including a Lindbladian description of finite-duration measurement 



C. A. Brasio and R. d. J. Napolitano 

Instituto de Fisica de Sao Carlos, Universidade de Sao Paulo, 
P.O. Box 369, 13560-970, Sao Carlos, SP, Brazil 

Abstract 

We consider the problem of the measurement of a system occurrmg during a finite time interval, 
while environmentally-induced noise decreases the system-state coherence. We assume a Markovian 
measuring device and, therefore, use a Lindbladian description for the measurement dynamics. For 
studying the case of noise produced by a non-Mar kovian environment, whose definition does not 
include the measuring apparatus, we use the Redfield approach to the interaction between system 
and environment. In the present hybrid theory, to trace out the environmental degrees of freedom, 
we introduce an analytic method based on superoperator algebra and Nakajima-Zwanzig projectors. 
The resulting master equation, describing the reduced system dynamics, is illustrated in the case 
of a qubit under phase noise during a finite-time measurement. 
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I. INTRODUCTION 

Quantum dynamical semigroups ^ constitute a powerful mathematical formalism for the 
treatment of non-unitary processes in quantum mechanics, which leads to many applications 
in the field of open systems, such as analysis of dissipation, decoherence and quantum 
measurement theory 2i]. It has been emphasized 5|, ISj that such phenomena can not have 
an adequate description by the conventional formalism of quantum mechanics, based solely 
on the Liouville-von Neumann equation. 

The most general form for the generator of quantum dynamical semigroups is the 
Lindblad[3, sl, with several applications to Markovian processes ^ and, more recently, 
to non-Markovian situations as well jo], |lo|. Indeed, as demonstrated in even in the case 
of the quantum Brownian motion it is possible to transform the Caldeira-Leggett master 
equation {ll| into a Lindblad equation, with the addition of a term that becomes small in 
the high-temperature limit. Moreover, in its stochastic form, the Lindblad e quat ion can be 
used in numerical methods, such as in the quantum-state diffusion approach 
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Thus, the Lindblad equation is widely used to describe irreversible quantum phenomena, 
such as the theoretical representation of finite-duration measurements. If we assume that the 
environment of an open quantum system under measurement does not include the measuring 
device, the measurement can, therefore, be approached using a Lindblad equation involving a 
set of Lindbladian operators taken as the measured observables, but with the environmental 
degrees of freedom also included in the global description. The next step to be taken in 
such a model is to trace the environment variables out of the formulation, leaving only 
a master equation describing both, the system and the set of measuring observables. The 
tracing procedure is usually very complicated due to the non-commutativity of several terms 
in the total Hamiltonian comprising the system under scrutiny, its environment, and the 
observables to be measured. 

Moreover, in the finite-time measurement, we can consider a Markovian interaction be- 
tween the system to be measured and the measuring apparatus. In that case, the relevant 
specification of the apparatus is effectively contained in the Lindbladian and the problem can 
be treated with the usual tools of that context. However, if we consider that the interaction 
between the system and its environment is non-Markovian, the action of the environmental 
degrees of freedon will be described by the Liouvillian only. 
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In this case, the entanglement between the relevant system and its environment compli- 
cates the description of the reduced time evolution of the system. One way of taking into ac- 
count the correlations between the system and its environment involves the use of a projector 
superoperator that transforms the total density operator p into its reduced counterpart, ps, 
calculated at the instant under consideration, tensorially multiplied by the reduced density 



where the most popular approach is the one by Nakajima and Zwanzig 
motivated by the work of van Hove, Prigogine, and Resibois 
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20] ■ Even so, however, it is 



still necessary to trace out the interaction to obtain the temporal equation for the relevant 
part of the system. Accomplish this task using directly the exponentials of each term to 
suppress them makes the calculations extremely difficult and prohibitively expensive. The 
procedure so hard, giving rise to many errors. 

In this paper, we will analyze finite-time measurent for a case of Markovian interaction 
between the principal system and the measurement and non-Markovian interaction with 
the system and the environment. We present a method to obtain the master equation in 
the Born-Markov approximation, by tracing out the environmental degrees of freedom using 
the time-independent thermodynamic projectors and superoperators for each Hamiltonian 
(system, environment and interaction), including the Lindbladian term. The method is 
easily demonstrated and the superoperator properties substantially simplify the calculations, 
rendering the resulting master equation compact and very simple. 

Although we use the Born-Markov approximation and time-independent thermodynamic 
projector superorerators, we believe it to be possible to generalize the present approach 



22j, non-Markovian cases 
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26|, or 



28|. 



to higher orders of the perturbation series, as in |2l[. 
time-dependent thermodynamic projectors, as inj27l. 

The paper is structured as follows: the general problem is formulated in the Sec. 2, with 
the definition of superoperators and thermodynamic projectors; in Sec. 3, the necessary 
properties of these superoperators are demonstrated; and, in Sec. 4, the master equation is 
obtained. Finally, in Se. 5, we illustrate the method in the case of a single qubit, subject to 
phase errors induced by the environment, when a finite-duration measurement is performed 
on an observable that commute with the unperturbed qubit Hamiltonian. 
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II. DEFINITIONS 



The most general form for a master equation, according to the quantum dynamical semi- 
groups approach, is the Lindblad equation {l, 2, 7, 8|, 



d ^ i 



H, Psb] + Yl [^^PsbL] - I {lILjJsb]] , (1) 

where H is the total Hamiltonian and the Lj are the Lindblads (Hermitian operators for mea- 
surement description, not-Hermitian for dissipation description) that represent a Markovian 
interaction. Here we will consider the study of finite-time measurement. 

It is important to point out that the Lindbladian operator of Eq. ([1]) appears as an 
effective device to emulate a measurement on the system. Thus, the Lindbladian evolution 
plays the role of the measuring-device action on the system. In fact, had we started from 
an ab initio microscopic description of a measuring apparatus, we would have to distinguish 
between its microscopic and macroscopic variables. For each of the possible eigenvalues of the 
observable being measured, there must correspond a single value of a macroscopic variable. 
However, such a variable alone cannot completely describe the state of the apparatus, for 
there is a large number of its microscopic states that correspond to the same value of the 
macroscopic variable. Since, in a measurement, the only variable that matters to the observer 
is the macroscopic one, we can trace out the microscopic degrees of freedom of the apparatus. 
This procedure results in a master equation which, under the Born-Markov approximation 
{2!, asssumes the Lindblad form, which we adopt here 0, li]. 

Let us consider that the Hamiltonian H can be written as 

a = Hb + HsB + Hs, 

where Hb describes the environment, Hs describes the relevant system, and Hsb is the 
Hamiltonian for the system-environment non-Markovian interaction. Here we assume, as 
usual, that the environment-system interaction is factorizable, i.e., 

Hsb = (2) 

k 

where, for each k, operates only on the system and B^, only on the environment. The 
form of the interaction, Eq. ([2]), is general, satisfied by both amplitude and phase damping 



models. By hypothesis, the Lindblads Lj will operate only on the system S. Our aim is to 
obtain an equation for the time evolution of the system reduced density matrix, ps, 



Ps (t) = tVB {pSB {t)} . 
Let us begin by defining some superoperators. For any density-matrix operator X, 



BX = - 



h L 



Hb. X 



(3) 



SX 



Hs^x]+Y.{LA^-\{L'ux]y 



(4) 



and 



FX 



h 



HsB, X 



(5) 



We will also use the Nakajima-Zwanzig thermodynamic projectors V and Q 
V is such that its action is defined, for any density operator pB (to) i by 



It is easy to check that 



We also define 



PXit) = pB{to)^TTB<X{t) 



V = V. 



(6) 



Q = x-r, 



(7) 



where X is the identity superoperator {XX (t) = X {t)). It follows from that definition that 
Q is also a projector, i.e., 

= Q. 

III. THE SUPEROPERATOR PROPERTIES OF B, S AND T 

As a consequence of the definitions (|3]) and (j4]), B acts only on the environment and S, 
only on the system. Hence, 



BS = SB. 
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It follows from that commutation relation that 

exp {St + Bt) = exp (St) exp (Bt) = exp (Bt) exp (St) . 

For simplicity, let us consider the initial time as zero (to = 0) and that the global density 
operator is initially factorized: 

PSB (0) = PS (0) ® PB (0) . 
For the sake of simplicity, let us ignore the ® symbol henceforth and write 

PSB (0) = PS (0) PB (0) = PB (0) PS (0) . 

The partial trace, over the environmental degrees of freedom, of the resulting action of 
exp {—Bt) on the global density operator, can be written as 

TTB{exp{-Bt)psBm = y2^-J-Tr^{B-psBm 

n=0 



Tr^ {psB (t)} + E S^Trs {B^psB (t)} • 



But, 



TrB{S>SB(t)} = TrB{i3i3""V5B(t)} 

Hb.B'^-^Psb (t) 



(8) 



The trace is basis independent and, therefore, it is convenient, for the treatment of Eq. ([8]) 
to use the environmental basis {|A;)} of eigenstates of Hb, i.e., 

HB\k) = Ei\k), 

for all k (where the eingenstats can be degenerate or not). With this consideration, 
Tib {B'' PSB m = "^E^^l [^b, ^"^^55 (t)] 1^) 



-- (A;| HbB^-'psb (t) \k) - {k\ B^-^psB {t) Hb \k) 
n 

~ {k\ E^B-'psB it) \k) - {k\ B-'psB (t) \k) 



[E^ B-'psB it) \k) - E^ {k\ B-'psB (t) \k)] 



0, for all k. 



Therefore, 



Tib {exp(-St) psB {t)] 



Ttb {psB m + y, ^^TrB {B-psB (t)} 

n=l 

TrsipsB (t)} + 



From this property, it follows that 



exp (St) ps (t) 



exp {-St) Ticb {psB (t)} 
= exp {-St) Ttb {exp {-Bt) psB {t)} 
= Ttb {exp {-St) exp {-Bt) psB {t)} 
= Ttb {exp {-St -Bt) PSB {t)}. 

For the interaction superoperator of Eq. ([5]), we have, from the ([2]) hypothesis. 



= - 



h 



SkBk, X 



for any density operator X. 



IV. TRACING OUT THE ENVIRONMENTAL DEGREES OF FREEDOM 



The equation that we want to solve is 



Hs + Hb + Hsb, Psb {t)] + Yl (^^JPsbL] - \ {^j^, P^b 



In terms of the superoperators of Eqs. ([3]) 



d 



and ([5]), the equation becomes 



-rJsB{t) = {s + B + :F)psB{t). 

at 

With the properties expressed by Eqs. ^ and fllOl) . let us define the operator 



a {t) = exp {—St — Bt) psB {t) ■ 



Then, 



ps{t) = exp{St)TTB{a{t)} 



and the equation satisfied by a (t) is written as 



d 
dt 



exp {-St - Bt) 



PSB it) + exp {-St - Bt) -j^PsB {t) 

d 



= -{S + B) exp {-St - Bt) PSB {t) + exp {-St - Bt) -psB {t) ■ 
We can insert Eq. (fTTi) into Eq. and the result is 

^a(t) = exp{-St-Bt)J^exp{St + Bt)a{t). 
Let us define, then, the superoperator: 



g {t) = exp {-St - Bt) J^exp {St + Bt) . 



(14) 



(15) 



We note that Q {t) is of the order of magnitude of the interaction Hsb, because it is propor- 
tional to the super-operator T that, in turn, is of the order of magnitude of Hsb- Thus, we 
have 

d 



dt 



a{t) = g{t)a{t) 



(16) 



Up to this moment, we have used superoperators based on each of the Hamiltonian and 
Lindbladian terms {B, S e T) and obtained Eq. (fT6|) . which involves the global density 
operator psB {t)- Now, let us use the Nakajima-Zwanzig thermodynamic projectors, defined 
by Eqs. §^ and in Eq. ([T^: 



and 



^ d ^ , . 



vg {t) a {t) 

Vg {t) {V+Q)a {t) 

vg {t) Va {t) + vg {t) Qa {t) 

vg {t) V^a {t) + vg {t) Q^a {t) 

Qg {t)a{t) 

Qg {t) {V+Q)a{t) 

Qg {t) Va {t) + Qg {t) Qa {t) 

Qg {t) V^a{t) + Qg {t) Q^a{t). 



Since those projectors are time-independent, it follows that 



I [Va {t)] = [VG {t) V] [Va (t)] + [VG (t) Q] [Qa (t)] , 
i [Qa (t)] = [QG (t) V] [Va (t)] + [QG (t) Q] [Qa (t)] . 
Let us formally integrate the second of the Eqs. ( ITTl) : 

Qa {t) = Qa (0) + T dt' [QG (f) V] [Va (t')] + T dt' [QG (f) Q] [Qa {t')] . 



(17) 



From the definition of a {t) , Eq. f|T2|) . we know that 

a (0) = psB (0) 

= Ps{0)pb{0) 



and, therefore, 



Hence, 



Qa (t) 



Qa(0) = iX-V)psiO)pB{0) 

= PS (0) PB (0) - Vps (0) PB (0) 
= PS (0) PB (0) - PB (0) TrB {ps (0) PB (0)} 
= p5 (0) h (0) - (0) Ps (0) Trs {pB (0)} 
= Ps{0)pb{0)-Pb{0)ps{0) 
= 0. 

dt' [QG (f) V] [Va it')] + / dt' [QG (f) Q] [Qa {t')] 



(18) 



'0 ^0 

showing that Qa (t) is of the order of magnitude of Hsb- Here we use the Born approximation 
and only keep terms up to the second order of Hsb- Accordingly, the second-order iteration 
of Eq. (fT8|) gives 



Qa (t) 



dt' [QG{t')V] [Va{t')] 



+ 



dt' [QG (f) Q] 



dt" [QG it") V] [Va it")] 



(19) 



Let us formally integrate the first of the Eqs. ( ITTI) : 

Va it) = Va (0) + [ dt' [VG it') V] [Va it')] + [ dt' [VG it') Q] [Qa it')] . (20) 
Jo Jo 
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Substituing Eq. (l20l) in the first of tlie Eqs. (fT7|) . we obtain 
d 



[Pa {t)] = [Vg (t) V] [Va (0)] + [VG (t) V] / dt' [VQ {t') V] [Va {t')] 
dt Jo 

+ \V<3 it) V] f dt' [Vg (f) Q] [Qa it')] + [Vg it) Q] [Qa it)] . (21) 
Jo 



From Eq. ([T9|) we see that the third term on the right-hand side of Eq. (12T|) is of the 
third order in Hsb, and, therefore, we neglect it. Substituing Eq. ( fT9l) in the forth term of 
right-hand side of Eq. f|2ip and keeping only contributions up to second order in Hsb, we 
obtain 

^ [Pa it)] = [pg it) p] [Pa (0)] + [pg it) p] dt' [pg if) p] [Pdit')] 

+ [Pgit)Q] f dt' [Qgit')P][Pait')] 
Jo 

= [pg it) p] [Pa (0)] + [ dt' [pg it) p] [pg it') p] [Pa it')] 

Jo 



'dt' [Pgit) Q] [Qgit')P] [Pa it')] 







[pg it) p] Pa (0) + / dt' [pg it) ppg it') p] [Pa it')] 

^ Jo 

dt' [Pgit)QQg it')P] [Pa it')] 

[pg it) p] [Pa (0)] + f dt' [pg it) pg it') p] [Pa it')] 



I.e., 



dt' [Pgit) Qgit')P][Pait')] 
= [Pg it) P] [Pa (0)] 

+ J dt' [pg it) pg it') p + pg it) Qg it') p] [Pa it')] 

= [Pg it) P] [Pa (0)] 

+ [ dt' [pg it) iP + Q)g it') p] [Pa it')] , 

Jo 

I [Pa it)] = [Pg it) P] [Pa (0)] + 1^ dt' [pg it) g it') P] [Pa it')] . 
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The first term of this equation can be written as 



[Vg (t) V] [Pa (0)] 



But, following from Eq. ([T5 



vg (t) VVa (0) 
vg (t) Va (0) 

vg (t) vps io)PB (0) 

vg (t) pb (0)TrB{ps (0)/)b (0)} 
vg it) PB (0)/)s (0)TrB{pB (0)} 
vg (t) Pb {0)ps (0) 
PB {0)TrB{git) PB (0) PS (0)} 
PB (0) Tib {git) PB mPs (0) . 



Trs {g (t) PB (0)} PS (0) = Tib {exp {-St - Bt) J^exp {St + (0)} ps (0) 

= TrB {exp {-St) exp (--Bt) J^exp exp {St) pb (0)} ps (0) 
= exp {-St) Trs {exp {-Bt) J^exp (Bt) (0)} exp {St) ps (0) , 



I.e., 



[Vg {t) V] [Va 



PB (0) exp {-St) TiB {exp {-Bt) ^exp (Bt) (0)} exp {St) ps (0) 



Now, let us use Eq. ([2]): 



TiB {exp (-Bt) ^exp {Bt) pb (0)} = Tr^ {exp {-Bt) [S^Bk, [exp (Bt) pb (0)]]} 



~ E {exp {-Bt) Sk [Bk, [exp (Bt) pb (0)]]} 



Analogously to the calculation leading to Eq. ([9]), we obtain 



^ 5^ 5fcTrB {exp (-Bt) [5^, [exp (Bt) pb (0)]]} 



TrB I exp (-Bt) Bk, [exp (Bt) p^ 



TrB I Bk, [exp (Bt) p^ 
TrB j-Bfc [exp (Bt) ps 
TrB|[exp (Bt) p^ (0)] 4 
0, 
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which imphes that 



[Vg (t) V] [Va (0)] = 

and 

j^[Va{t)] = f^dt' [VG{t)G{t')V][Va{t')] 

= f dt' [vg{t)g{t')vva{t')] 

Jo 

= [ dt' [Vg{t)g{t')Va{t')]. (22) 
Integrating Eq. (122|) between t' and t yields 

Va (t) - Va {t') = f dt" f dt'" [Vg (t) g {t'") Va {t'")] , 
Jt' Jo 

which shows that the difference between Va (t) and Va (t') is of the second order of mag- 
nitude in HsB and, therefore, we can write Va (t) instead of Va {t') in the integrand of Eq. 
( !22|) . obtaining an equation that obeys the Markov approximation, without violating the 
Born approximation. Thus, the master equation that we finally obtain, is written 

^ [Va (t)] = dt' [Vg (t) g (f) Va (t)] . (23) 

In the Born approximation, to obtain the reduced density operator of Eq. (IT^ . we must 
solve Eq. f l23ll using the initial condition 

TTB{a{0)} = TrB{ps{0)pBm 

= ps{0)TlB{pBm 

= PsiO). 

V. EXAMPLE 

As a simple example of our method, let us consider a two-level system (a spin | par- 
ticle pointing along the z direction) in contact with a thermal bath of quantum harmonic 
oscillators: 

Hs = hwoffz 
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and 



We take the interaction to be the phase damping where, in reference to Eq. 

Sk = ha^, 

^Bk =gkbi + glbk, 
giving 



HsB = h^a, {^Qkbl + glb,^ . 



(24) 



For the initial state of the thermal bath, let us consider the vacuum state (T = 0): 



PS = (|0)|0)...)®((0|(0|... 



(25) 



Now we study the case of a measurement of the z component, using a single Lindblad: 



L = Xa^ 



where A is a real number. Then, Eq. (jlj) gives 



SX = —iujQ 



}? ( CTzXcTz — X 



To simplify the notation, let us define the following quantities: 



R (t) = exp i-St) PS (t) 



(26) 



and 



Va (t) = R{t) pb. 

The action of exp [St) and exp{Bt) can be calculated in the following way. For an arbitrary 
density operator ^(0), let us define 
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X (t) = exp {St) X (0) . 



Hence, 



di 



X{t)=S exp (St) X (0) = SX (t) 



that is, 



d 
di 



X (t) = —iuo 



a.Xa. - X 



The solution of Eq. (128|) can be easily determinated 13| : 



(27) 



(28) 



(29) 



Xn{t) =Xn(0), 

Xi2 (t) = Xi2 (0) [cos {2uot) - isen {2uot)] . 

Analogously to the case of exp (St), for an arbitrary density operator X{0), let us define 



X{t) = exp {Bt)X{0). 



From Eq. ([3]), it follows that 



HB,Xit) 



whose solution is, simply, 



X (t) = e-'^'X (0) e'^\ (30) 
From the above explicit actions of exp (St) and exp(i3t) it easily follows that 



exp (St) exp (Sf) = exp [S {t + t')] , 

exp (Bt) exp (Bf) = exp [B {t + t')] . 

With Eqs. (130|) and (129|) we are able to solve Eq. ( 123|) . Separating the system and 
environment terms, we obtain, in terms of Eq. fl26|) : 
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Vg {t)G{t')Va it) = 



PB + 



_g-5t |g<S(t-t') 



E4' (e^*'pB) 



. k 

Expanding the environmental superoperators using Eq. (l30l) and grouping the similar terms, 
we have: 



vg {t)g{t')Va (t) 



e^^'R it)) a. 



St I „S{t^t') 



xtrs <j ft *E^fce 



+ e 



-St } ^S(t~t') 



St; 



k' 

oS{t-t') 



(31) 



e 



pB X 



k k' 



h 



Using the initial state of the environment, Eq. fl2^ . and the interaction between the system 
and the environment, Eq. (Elj), we trace out the environmental degrees of freedom and 
write: 



vg {t)g{t')va (t) 



{e~^'a, je'^C*-*') [(^e'^^'i? (t)) a,] | - e^-^* |e<5(*-*') ^(^e^*' R (t)) a,] } a,} ® p^? x 



xE \9kf {cos [uk {t - t')] + isen [u^ {t - t')]} 



+ 



,S{t~t') 



aJe^''R{t) 



xE \9k\^ {cos [uk {t - t')] - isen [uk (t - t')]} 



pB X 



Now we rewrite the system superoperators using Eq. ( l29i) . Thus, in terms of R(t) , we 
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have 



vg {t)g{t')va it) = 

/?12 \ 2 

E \9k\ COS [uk {t - t')] , 
R21 k 



which, in accordance with Eq. (l23|) . gives 



d I Ru R 



12 



R21 R22 J \ R21 



4 1 ° ^''] f dt'y2\g,fcos[u,{t-t')]. (32) 

\ Ro^ Jo V 



Now we make the continuum transformation Let us define the density of states as 

k 

which allows us to rewrite Eq. (132|) as 



d I Ru R12 \ . I ^ Ri 

di 



R21 R22 I \R2i I Jo Jo 

and, with the change of variable t = t — t', we obtain 



R12 \ , r°° 

4 \ dt' duJ{uj)cos[uit-t')] 

\ /?oi n / Jo Jo 



d I R^^ ^12 I I ^ R\2 . , 

' ' 4 \ duJ{u)cos{uT). (33) 



dt \ i?2i i?22 ) \R2i J Jo Jo 

To obtain the solution for the diagonal terms of R (t), we do not need further considera- 
tion: 

Ru{t)= Ru{0) 

R22{t)= R22{0) 

However, for the non-diagonal terms, it is necessary to specify J {u) . Let us use the Ohmic 
density of states: 

J (u) = riuje~'^ , 

where r] and Q are real and positive constants. Hence, Eq. (155]) becomes 



d I Ru R 

dt 



12 



R21 R22 j \ R21 
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( Ri2\ /■* /""^ „^ 
4:7] \ \ dr dojcoe " cos (ut) 

\ i?9i / ^0 Jo 



The double integral gives 



dr / dujue^'^cos (oor) = - 



^0 



I.e., 



dt 



dt 



R 



12 



R 



21 



The solution for this system is 




Rl2(0) 
fl2l(0) 



At last, now we calculate the density-operator elements by inverting Eq. ( l26l) and we obtain 
Pii{t)= Pii(O), 



(34) 

- isen [2ujQt)\ , 

remembering that P22 if) = ^ — Pii f) and p2i (0 = P12 (0 • 



P12 (t) = P12 (0) ^^^^^^^2^^ [cos {2uot) - isen {2ujot)] 



VI. CONCLUSION 



In summary, here we present a new method to describe the dynamics of measurements 
that occur during a finite time interval, while the system being measured interacts with the 
rest of the universe and, due to the consequent environment ally- induced noise, undergoes de- 
coherence. We use a Lindbladian description of the measuring apparatus, whose interaction 
with the system we assume as Markovian. To treat the noise introduced by the fact that, 
during the finite-duration measurement, the system is perturbed by the environment, we use 
a Redfield approach to the dynamical description of the interaction between the system and 
its environment, assumed non-Markovian. The resulting unprecedented hybrid description 
is shown to be capable of substantially simplifying the tracing procedure, which is usually 
very complicated due to the non-commutativity of several terms in the total Hamiltonian, 
comprising the system under scrutiny, its environment, and the observables to be measured. 
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The superoperators defined in Sec. 2 introduce simplifications of the calculations leading 
to Eq. fl23|) . that is compact and can be solved in terms of the unperturbed solutions (in Sec. 
5, for example, we used Eq. (129|) ). Moreover, for the Born-Markov approximation, regardless 
of the model chosen for the environment, the reduction of the density operator becomes 
evident, as can be verified in Eq. (13 ip . The simple phase-damping-interaction example 
of Sec. 5 (see Eq. (El])) at zero temperature already provides an important and expected 
result in the quantum mechanics of open systems: the intensification of the environmentally- 
induced decoherence, as indicated in the denominator of the second of Eqs. (134|) . evidencing 
the power and convenience of the present approach. 
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